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We study the non-singlet sectors of Matrix Quantum Mechanics in application to two-dimensional 
string theory. We use the chiral formalism, which operates directly with incoming and outgoing 
asymptotic states, related by a scattering operator. We argue that a general non-singlet scattering 
amplitude decomposes into scattering amplitudes in the adjoint. In the adjoint representation, 
we express the phase of the tree-level scattering amplitude as an integral over the Fermi sea. We 
claim that this expression holds for any, in general time-dependent, profile of the Fermi sea. In 
the case of stationary Fermi sea our formula reproduces Maldacena's scattering amplitude for a 
long string to go in and come back to infinity. 



1. Introduction 

The two-dimensional string theory arises as the collective field theory of the gauged matrix 
quantum mechanics in "upside-down" gaussian potential, or shortly MQM [|1].0 MQM 
involves one gauge field A\ and one scalar field X?, both hermitian N x N matrices. The 
theory is formally described by the action 



S = J dt Tr [PV A X - \{P 2 -X 2 )] , 



(1.1) 



where VaX = dtX — i[A,X] is the covariant time derivative. Since the potential is 
bottomless, one should introduce a cutoff and tune the size iV with the cutoff before 
taking the large iV limit. 

The tachyon dynamics in 2D string theory with flat, linear dilaton, background is 
described by the singlet sector of the Hilbert space. In the singlet sector the gauge field 
simply plays the role of Lagrange multiplier assuring that the matrix coordinate and mo- 
mentum can be simultaneously diagonalized. As a consequence, the action (|1.1| ) describes 
a system of N non-relativistic free fermions in the upside-down quadratic potential. The 
ground state of the system is characterized by the Fermi level Ep = —ft, where g s = 1/fj, 
is the string coupling constant. 

The states obtained by introducing one or more Wilson lines in the adjoint represen- 
tation belong to the non-singlet sector of the Hilbert space. The wave functions in this 
sector transform by non-trivial representations of SU(N) that can be obtained by ten- 
soring the adjoint. Besides tachyons the non-singlet sector contains excitations which in 
the compactified Euclidean theory correspond to winding modes, or vortices, on the world 
sheet HHjjU- The non-singlet sector may contain more general string backgrounds. It was 



conjectured |J [7] that adding two oppositely oriented Polyakov loops in the action fll.l 




SS = ATrexp ( / dtA(t) J + ATrexp I - dtA(t) ) , (1.2) 

may deform the metric of the target space and produce Euclidean black hole background 
||. It is known that the deformations by winding modes are integrable |J and the de- 
formed partition function is a tau- function of a non-compact Toda hierarchy |7]] []1(J . How- 



ever, knowing the partition function is not sufficient to reproduce the geometry of the 



1 The relation between the 2D string theory and MQM was originally proposed by Kazakov 
and Migdal §. 



1 



background. In order to probe the background curvature, one should be able to calculate 
scattering amplitudes of tachyons in presence of a deformation by a non-singlet source. 
This problem is well posed in the Lorentzian formulation of the theory. 

Recently, Maldacena gave a world-sheet interpretation of the non-singlets in 
Lorentzian MQM. He considered mainly the states containing one Wilson line, which 
transform according to the simplest non-trivial representation, the adjoint. According 
to Maldacena, the wave function in the adjoint describes a particle-like excitation (impu- 
rity), interacting with the Fermi sea. The collective coordinate of this "adjoint particle", 
describes the position of the tip of a folded string that stretches from infinity. This in- 
terpretation passed several consistency tests. The tree level scattering amplitude for the 
tip of the long string, evaluated in |11| using the FZZ boundary Liouville amplitude fl2f 



turned out to be the same as the one extracted from the wave function in the adjoint 
representation, TZij.^i(Q) = fi^C^ — o~ik&ji subsequently calculated in [O. Furthermore, 



the density of states corresponding to the phase shift reproduced correctly the vortex- 
antivortex correlation function in the compactified Euclidean theory Further study 
of long strings in 2D string theory is presented in |14| . 



In order to make the next steps towards constructing a solvable model of Lorentzian 
2D black hole it is important to extend this picture to the whole non-singlet sector. As all 
allowed representations are contained in the tensor products of the adjoint, all non-singlet 
states can be described in terms of tachyons and adjoint particles, or long strings. The 
wave function of such a state can be in principle computed by solving the corresponding 
matrix Calogero equation. However, in spite of the fact that problem is integrable ||16|| , 
this seems to be a rather hopeless task. Some interesting speculations concerning higher 
representations were presented in ||17|1 , but no quantitative description exists so far. 



In this paper we argue that the scattering problem in the non-singlet sector of MQM 



can be solved by means of the chiral formalism introduced in |18| and later developed in 
I^pD|[^P|P^1[rr4|[p4]|2ip[p5[P^, in which one performs a canonical transformation to 



the light cone variables 



X 



+ 



x+p 

V2 ■ 



X. 



x-p 

V2 ■ 



(1.3) 



The wave functions in X + and X spaces describe the asymptotic states in the original 
theory respectively in the infinite past and future. The outgoing and the incoming states 
are related by Fourier transformation 



;i-4) 



which represents the scattering operator in the chiral basis. The Hamiltonian correspond- 
ing to the action ( p. . 1|) is bilinear in the new variables, 



H = -\ r Ti(X + X_ +X_X + ), (1.5) 
and the general solution of the corresponding Schrodinger equation is 

$ ± (X ± ,t) = e^ 2 * $±( e ^X ± ). (1.6) 

In the non-singlet sector the wave function is a vector transforming according to a 
(not necessarily irreducible) representation 1Z of SU(N): 

dim7?. 

$f(nx ± rt) = $J(xj, n e su(N), (1.7) 

j=i 

where 1Zi^j{Q) is the matrix of the representation. The U(N) symmetry allows to reduce, 
as in the singlet sector, the original iV 2 degrees of freedom to the N eigenvalues xf . . . x^- 
of the matrix X + or X . 

Our claim is that one can construct a large iV collective field theory for the scattering 
amplitude between the states $~ and $ + , which is given by the inner product 

/dim7?. 
dX + dX_ e lTrX + x - ®i( X -) ®t( X + )- (1-8) 

j=i 

This is achieved in two steps. The first step consists in integrating out the angular degrees 
of freedom in the matrix integration measure 

dX ± = dtt± dxf... dx% A 2 (X ± ), (1.9) 

where A(X ± ) = Yli<j( x t ~ x f) * s the Vandermonde determinant. The second step is 
to formulate the result in terms of the collective field, the phase space eigenvalue density 
p(x + , x~). 

As a consequence of ( |1.7| ), the angular dependence of the integrand in ( |1.8| ) is only 



through the factor 

dim7^ 
K=l 



The integral with respect to O = fi+fil depends only on the representation 1Z and not on 
the wave functions: 



i?Ax +1 x_) 



ix , nx n f 



;i.io) 



SU(N) 



In the singlet sector, 1Z = I, this is Harish-Chandra-Itzykson-Zuber integral [p8|. The 
general integral (|1.10|) has been first studied by Shatashvili |29|] who found the complete 
solution, but in a form not explicitly invariant under permutations of eigenvalues and hence 
not immediately applicable in the large N limit. A nice symmetric formula for the adjoint 
representation, lZikji(£l) = Qij£V lk — fiifcOjj was originally proposed by Morozov |3(| and 



subsequently proved by Eynard and collaborators [j31j|3^] . Generalization of Morozov- 
Eynard formula for any representation was found in Therefore the first step, the 

U(N) integration, is essentially accomplished. 

In order to extract the large iV limit and construct the collective field theory, we 
need to evaluate the inner product ( |1.8| ) for the wave functions that describe low energy 
excitations above the ground state. The latter is constructed by to filling the fermionic 
energy levels of the singlet sector up to the the Fermi level Ep = It is sufficient to 
solve the problem for a standard set of functions spanning the sector containing p adjoint 
particles, 



$± . . . 

*1 " "*p iJl " 'Jp 



(£l J "•! ^p'l 



1 




1 




l l-7(T±(l) 





d>±(X ± ), (1.11) 



CT =t( P ) 



where $ is the ground state wave function and a ± are permutations from Sn- Here we 



suppressed the time dependence, which can be reconstructed from ( |1.6|) . The functions 



( |1 . 1 1|) transform according to the tensor product of p adjoints. The projection to the 
irreducible components of the tensor product is done by summing over the permutations 
with a character of the symmetric group. 

The integral for the inner product of functions of the type ( |1 . 1 1|) is similar to that for 
the mixed correlators in the two-matrix model, which were evaluated in the large iV limit 
by Eynard and Orantin |34]]. The authors of [34| found that a generic mixed correlator 
decomposes as a sum of products of two-point mixed correlators. An important for us 
fact is that the coefficients in the sum do not depend on the radial part of the matrix 



measure. Therefore the result of 34 must be applicable also in the case of the inverse 
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matrix oscillator, which means that any scattering process involving only adjoint particles 
can be decomposed into one-particle scattering amplitudes. 

In this way the success of the chiral approach is guaranteed if it works for the adjoint 
representation. Our aim here is to demonstrate that this is indeed the case. We will obtain 
a formula for the scattering amplitude of the adjoint particle in terms of an integral over 
the Fermi sea, which we expect to be valid for any, in general time-dependent, tachyon 
background. In the case of stationary Fermi sea we recover Maldacena's expression |Ll] 
for the scattering amplitude for a long string to come in and go back to infinity. Our 
more general answer can be used to calculate the tree-level processes involving one adjoint 
particle and any number of tachyons. Scattering amplitudes involving several adjoint 
particles will be considered in a elsewhere. 



2. The scattering amplitude in the adjoint representation 



It will be convenient to absorb into the wave function a Vandermonde determinant from 
the measure (|1.9|): 

*±(X ± ):=A(X ± )$±(X ± ). (2.1) 
Then a complete set of wave functions in the adjoint sector is given b 

1 



ij 



det 

ki 



(2.2) 



The last factor is the general eigenfunction of the matrix Hamiltonian in the singlet sector, 
which is a Slater determinant of one-fermion eigenfunctions 



tp%{x ± ) 



/2ir 



(x ± ) ±lE - 



(2.3) 



We are going to consider only states which represent incoming leftmovers and outcoming 
rightmovers. Then the eigenvalues of X ± can be assumed positive and the operator func- 
tions (|2.2|) are analytic in the complex £-plane cut along the negative real axis. One can 
think of £ + and £~ as the phase space coordinates of the adjoint particle. The Hamiltonian 
( pT5|) acts on the wave function (|2.2j ) as 



± d 



+ l) + Ef + ... + E 



(2.4) 



2 Since we do not subtract the trace, this function has a small component in the singlet sector. 
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We will be eventually interested in evaluating the inner product for the eigenfunctions of 
the Hamiltonian, 



#! ± (XJ = (XJ±^ det ^% ± (xf) , (2.5) 



A;2 



,± („±^ 



which are obtained as Mellin transforms of ( |2.2| ): 



^ ± (x ± ) = ±^^ rdU^ E± ^;x ± ). (2.6) 



o 



Our aim is to evaluate the scattering amplitude for the adjoint particle in the large 
N limit, when the eigenvalues can be replaced by a collective mean field. For that it is 
sufficient to consider the diagonal elements of the scattering operator, E^ = Ek, k = 
1, . . . , N, which implies E + = E~ . The non-diagonal elements, which are of subleading 
order, allow to study the back reaction of the adjoint particle to the collective field. 

We will proceed as follows. First we perform the angular integration in the inner 
product of the functions (|2.2| ) using Morozov-Eynard formula. Then we consider the large 
N limit, in which the result takes the form of Fredholm determinant, and evaluate the 
leading contribution at large cosmological constant \i. Finally we perform the Mellin 



transform according to (|2.6| ) and obtain the reflection factor as a function of the energy E. 



2.1. Eliminating the angles 

The angular integration in the inner product of the wave functions (|2.2|) is evaluated by 
Morozov-Eynard formula |3^] 

detf S+ij } v 1-detS 



V?_+x_ 12 e + +^ + y ~~ iA(x_)A(x + ) > 



SU(N) 



where S is the matrix with matrix elements Sjk = e lx i Xk and X ± = diag^j , . . ■%%)• 
Plugging Q2.7| ) in the definition of the inner product (|Q|), we get 



POO 

= -i / d N x+d N x- det (*,")] det [^ h (x+ 

J 



x det 



S + i , } v S 



(2A 



Here we dropped the second term in the Morozov-Eynard formula, which is needed to cancel 
the constant term in the expansion of the r.h.s. at £± — > oo. This term does not depend on 
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£ ± and therefore is irrelevant for the scattering phase. Note that the Vandermonde factors 
from the measure disappear due to the redefinition (|2.1|) of the wave function. 

A form of the inner product suitable for taking the large iV limit is obtained if we 
rewrite (|2.8|) as a determinant of double integrals: 

N 

(^-\^+) = R^(( + ,U ]jR(E k ), ( 2 -9) 

k=i 

R ad (H + iU = -idetiSiEj-E^+iKiE^Ek)} , ( 2 .10) 

jk 

K(E',E)= ( E'\ WT ^-^- ) \E). (2.11) 

Here R(E) denotes the fermion reflection coefficient, or the bounce factor, determined by 
the inner product of the one-particle eigenfunctions (|2.3|) , 



/•OO 

/ dx + dx_^ E J (x~) e lx+x -^' (x + ) =R(E) d(E-E'), (2.12) 
Jo 

and the r.h..s. of ( |2.11|) is defined as 



(E'lflE) — ^ dx + dx -^-'( x )f(x_,x + )e ix + x -^'(x+) 

\/R(E) R(E') 

In particular, 

(E'\1\E) =8{E'-E). (2.14) 



The second factor in the inner product (|2.9|) is the inner product of the singlet states whose 
wave functions are given by the Slater determinants in ( |2.2| ). Therefore the first factor, 
-R ad (£ +5 £_) 5 gives the scattering amplitude (more precisely, the reflection coefficient) for 
the adjoint particle as a function of the collective coordinates, £ + and , of the initial and 
final states. 

2.2. The scattering amplitude as a Fredholm determinant 



The expression ( j2.10| ) for the inner product was derived for a generic wave function of the 



form ( |2.2|) . Now we will specify the Slater determinant in (|2.2| ) to be the wave function of 
the Fermi sea filled up to Ep = —fJ>- For this purpose we introduce a large cutoff A ^> [i, 
say, by putting a wall at x ± ~ vA~, so that the spectrum becomes discrete. Then the 
density of the energy levels in the singlet sector is given by jl8| [|19 



^=^"^ UE) = -Ho 6 R(E), (2.15) 



and the level spacing is 



AE 



2n 



1 



(2.16) 



p(E) log A 

The total number N of eigenvalues is given by the number of energy levels between the 
bottom of the regularized Fermi sea, -Ecutoff = —A, and the Fermi level Ep = — fi: 

dEp(E). (2.17) 



N = 

Finally, the inner product (|2.12|) is replaced by the following relation for the regularized 
eigenfunctions: 

/•OO 

/ dx + dx_^ B J e lx + x -^l k = R(Ej) p(Ej) S jk . (2.18) 
Jo 

We would like to evaluate the scattering amplitude ( |2.10| ) for the case where 
Ei, ... , En are the allowed energy levels in the interval [—A, —fi\. In presence of a cutoff, 
the expression ( |2.10| ) takes the form 

# ad (£+,U = -idet [piE^Sjk + iKiE^Ek)] 

(2.19) 



jk 
N 



= -i TT P( E k) x det 5 ik + K(Ej, E k ) 

jk 



k=l 



The first factor is an irrelevant infinite constant and will be neglected. The second factor 
has a smooth limit A — > oo: 



irx,a = i+*£^A'(£^)-i£4^ 

3 j,k 



K{Ej, Ej) K(Ej, E k ) 
K(E k , Ej) K(E k ,E k ) 



+ 



dE 



1 

2! 



-/■« 




/ dE 


/ dE' 


J 2W 


J 2W 


-oo 


— oo 



K(E,E) K(E, E') 
K(E',E) K(E',E') 



+ ■ 



(2.20) 

The series is by definition the Fredholm determinant of the kernel ( |2.11| ) , restricted to the 
interval [— oo, — fj]: 

R ad (Z + ,U = Det[l + iK}. (2.21) 



Therefore the scattering phase S (£, ,£_), defined by 

i? ad (^ + ^_)=exp [iS(Z + ,U] 
has the following integral representation: 

s(U,L) = J2- J — TrKn ' 



n=l 



n 



[2.22) 



(2.23) 



TrK n. = I dE± . _ _ dE^_ x{E\, E 2 ) ■ ■ ■ K(E n , Ei). (2.24) 

J — oo 

2. 3. Tree level calculation of the reflection factor 

The Fredholm kernel K(E, E') can be evaluated semiclassically for large negative energies. 
The easy calculation is given in Appendix A. The result is 



K(E + e, E - e) 



V e { +E ) ( +E ) (2.25) 



i_J zsinh(27re) LL + E 



The corrections are of order 1/E. It is obvious that the integral for Tri^ n is convergent 
for n > 1 and behaves as o(^u 1_n ). Therefore at tree level the phase ( |2.23|) is given by the 
first term of the series, 

S((+,Z-)=TrK= [ ^ dE K(E, E), (2.26) 

J —A 

K(E,E)= ^fl-^- E \ (2.27) 

This integral is logarithmically divergent at minus infinity. To make it convergent, we cut 
it off at E = —A, where /i < A C A. The new cutoff A will be given below a precise 
meaning in terms of the world-sheet theory. Let us introduce the variable 

8= ilog(£ + £_/V) (2.28) 



and the function (it already appeared in Maldacena's paper) 

tv J ^ (tanh£ 



f(x) = - [ dt(-^— + (). (2.29) 



- oo 



Then we can write the regularized integral (|2.26| ) as 



S(£+, £-) = /(§ MA/V) - a) - / (s) 
^i(ilog(A/^)- S ) 2 -/(- S ) 



(2.30) 



In the second line we used the asymptotics f(x) ~ x 2 /tt at x 3> 1. 

The leading corrections to this tree-level formula come from the terms of order 1/fj, in 
the semi-classical expression for the kernel ( |2.27| ), as well as from the second term, TrK 2 , 
of the series ( |2.23| ) . It however natural to expect that the \j \x corrections cancel and the 



perturbation series for the scattering phase is in 1/fi . 
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2-4- The scattering phase as an integral over the Fermi sea 

In the derivation of the scattering phase we did not use the particular form ( |2.3|) the one- 
particle functions. Therefore the derivation remains valid also in presence of a tachyon 
source, whose only effect is that the one-particle wave functions get deformed [BJ. In this 
case the scattering phase will be given by the same expression (|2.23|) , but with a deformed 
kernel K. The leading term, the trace (|2.26|) , can be expressed alternatively as an integral 
over the Fermi sea, 

oo 

f dx.dx p(x,, x ) 

o 

where p(x +: x_) is the semi-classical density of the fermionic liquid. The expression ( |2.31| ) 
can be useful if we want to calculate the scattering phase in more general, time-dependent, 
tachyon backgrounds. In the case of a stationary Fermi sea the density is given by 

p(x + ,x_) = { 1 if><V_<A; (2.32) 
1 otherwise. 

2.5. Evaluation of the scattering phase in the energy space 



The scattering amplitude for the eigenstates with given energy (|2.5[) is obtained by applying 
the integral transformation (|2.6|) to both arguments of -R ad (£ + ,£_): 



+ • 

oo 



i _ _ r°° 

R ad (E+,E-) = — sinh(7rE + )sinh(7rE_) / e iS G+'*-\ (2.33) 

For a stationary Fermi sea the r.h.s. depends on £ + and £_ only through the variable s 
defined in ( |2.28|) , hence the l.h.s. contains a delta function imposing the energy conserva- 
tion: 



i? ad (£ + , E_) = e~ ld 5(E+ - E 

/■OO 

e -i8^{E) _ 1 e 2nE iE+1 / dg e 2s(iE+l)+i± ( s _ 1 log £)* -if(-s) _ 



"• ! < ! 2-kE »r-4-i / ' . •2,,o>:+h+iii.s-ii. i; . fi — s) (2.34) 

In the last equation we replaced the sine function in front of the integral by exponent, 
which is justified for large positive energies, E ~ ^- l°g(A//x). 

We subtract, as in |]Tl|j , the energy gap that separates the singlet and the adjoint 
sectors and introduce the shifted energy! 

e = E + — Int. (2.35) 

2tt A v ! 



3 Note that our energy variable E differs from the variable e used in |fj"l|| by a cutoff-dependent 
constant, E = e + ^- log A. 
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Then we write the integral (|2.34 ) in the form 



/oo 
(fae 2,(fi+l)+i/(.)-i 1 r/6 ) (236) 
-oo 

2 

where we used the property f(—x) + f(x) = + f • The integral can be taken exactly 
using the remarkable fact (see Appendix B) that the function e 2s+l ^ s ^ reproduces itself 
after Fourier transformation. Thus we obtain for the phase factor 

e -i5»*(E) ^ I iE e ±(£ log ±) 2 e -if{-ne)-in/6 

(2-37) 

= A e i7r (2¥ lo s A ) x e - i7r ( £ -^F lo g A ) i7r X e i/(7re ^. 

We get for the scattering phase, neglecting the large cutoff-dependent constant, 

5 ad (E) = n(E-± log A) 2 - f(ne). (2.38) 

The answer is in accord with Maldacena's calculation of the scattering phase based on 
the FZZ formula for the boundary two-point function in Liouville theory, eq. (A. 4) in 



Appendix A of plj . Comparing the quadratic cutoff-dependent factors, we see that the 
cutoff A in the integral over energies should be identified with the square of the large 
boundary cosmological constant of the FZZT brane, A = ii 2 B . 



3. Discussion 

In this paper, we explained how to use the chiral formalism of MQM in order to evaluate the 
scattering amplitudes in the non-singlet sector. The scattering operator is given by matrix 
Fourier transformation and its matrix elements can be in principle evaluated by applying 
recently discovered generalizations of the Harish-Chandra-Itzykson-Zuber formula. We 
considered in detail the scattering process in which the incoming and the outgoing states 
are in the adjoint representation. The wave function in the adjoint describes a particle-like 
excitation interacting with the Fermi sea, which was identified in [II] as a long folded 
string stretching from infinity. 

We showed that the expression for the scattering amplitude such an 'adjoint particle' 
can be written in the form of a Fredholm determinant. The Fredholm kernel depends 
on the profile of the Fermi sea of the singlet sector. In the case of stationary Fermi 
surface, x + x~ = fx, we reproduced Maldacena's phase shift for a single long string fll 



Our tree-level formula ( [2.31| ) for the scattering phase as an integral over the Fermi sea 
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is actually valid also for time-dependent perturbations of the Fermi surface, e.g., by a 
tachyon source. This more general expression can be used to evaluate the amplitudes of 
emission or absorption of tachyons by the long string. 

The scattering amplitudes involving p adjoint particles are given by the inner products 



of the wave functions (|1 . 1 1 ) . In the case when (cr + ) 1 a = 1, the tree-level amplitude 



factorizes into p one-particle amplitudes. In the general case, when (cr + ) 1 a is a per- 
mutation with n cycles, the scattering amplitude factorizes into a product of one-cycle 



amplitudes. For each such amplitude one can use the result of |34[] to express it as a sum 
of products of one-particle amplitudes: 



(3-1) 



<xes fc 



i=l 



In this formula all the dependence on the eigenvalue distribution is contained in the one- 
particle amplitudes, while the coefficients are universal and related only to the U(N) 
integration. For example, the one-cycle-of-length-two amplitude is given by 



^(fi) £2"; £1 > £2 



tit - ami - & 



(3.2) 



One can evaluate in this way the most general scattering process involving any number 
of tachyons and adjoint particles. The corresponding asymptotic states are obtained by 



replacing in ( p. . 1 1|) the last factor, the ground state singlet wave function, with an excited 
singlet state. 

An alternative approach to study the scattering in the non-singlet sector consists in 
diagonalizing the scattering matrix given by the inner product ( |1.8|) . This is possible 
because the problem is integrable. In the case of the "upside-up" matrix oscillator such 
an approach, based on a spectrum-generating algebra that generalizes of the singlet 



sector, was developed recently by Y. Hatsuda and Y. Matsuo |35(| . In the "upside-down" 
case the generators of this algebra create non-singlet discrete statesi and can be used to 
write down Ward identities for the S'-matrix elements. 
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Appendix A. Quasiclassical calculation of the integration kernel 

Consider first the simplest problem, the calculation of the diagonal E' = E of the kernel 
( [2 .11|) for large negative energy, —E ^> 1. Introduce the parametrization 

x± = ^ e ±T , £ ± = ^pe ±CT (A.l) 

and evaluate the r-integral in the numerator of ( |2.13|) : 



ix , x 



e 



r°° dr 1 

dT / rUL r ^ e ir " (A.2) 

Jo V r r + P + 2 v /rpcosn(r - a) 

dr rl E e ir lnr-lnp 



'r r — p 

The remaining integral in 



00 dr „iE„ir In r — In p 



r°° dr_ 

K(E, E) = h fj^ (A.3) 



can be evaluated semiclassically. Up to o(-4) correction the integrals in the numerator and 
in the denominator are saturated by the same saddle point, r = —E. Therefore in the 
leading order 

K(E, E) = = ■"(-*)- "K+H (A . 4) 

-E - p -E - £ . £ 

The off-diagonal elements K(E, E') are evaluated similarly. The r-integral in this case 
gives 

'dr . = J.Ml|^^ (A.5) 

r + p + 2wrp cosh(r — a) r — p zsmh27re 



and the saddle point approximation of the integral in r yields (|2.25| ) . 
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Appendix B. Properties of the function ( |2.29|) 

1) Symmetry: 



/(-*) + f(x) = — + I . (B.l) 

7T 



2) Series expansion at x — > +00: 



/(*) = ^ + -- £(- + 4) e_2nx - ( B - 2 ) 

n>l 



3) Functional identity: 

(g-iTrS,, + e 2x^ e if(x) = e if(x)^ (B.3) 

4) Fourier transform: 

/oo 

ds e 2s e 2ls ^ +lf{s) = ie~ 2 ^ e - i/(_7rg) . (B.4) 
-00 



The symmetry property ( |B . 1| ) and the series expansion ( |B.2| ) follow directly from the 
definition ( |2.29| ). The functional equation is equivalent to 



f(x-in) - f(x) = -i dx(coth(x) + 1) = -zlog(l-e^). (B.5) 

J —00 

We do not have a complete analytic proof for the fourier transform ([B .4|) , so we had to 
complete our analytic argument with numerical evaluation. Equation ( [B.3] ) is symmetric 
under Fourier transformation 



— [ dye-- ipx+lf(y \ ehv) = I dse« ipx+if(x \ 



Therefore the Fourier image satisfies the same functional equation. There are two obvious 
solutions, f(p) = f(p) and f(p) = —f(—p), up to a periodic function under p p + i%. 
The saddle point evaluation of the integral, valid for large negative p, is compatible with 
the second choice, 

f(p) = -f(-p) -Hnh(p), 
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where h is periodic, h(p + in) = h(p). Further, the function h(p) must vanish exponentially 

2 

at p — > — oo, where the quasiclassics holds and we know that f(p) ~ — — . Therefore h is 
expanded at p — > — oo as a power series in e 2p . Finally, since f(x) vanishes exponentially 
for large negative x, the function h(p) must have a simple pole at p = 0. These 
condition fix the form of function h as 

. 1 + h ie 2 P + h 2 e 4 P + ■ ■ ■ 
hyp) = —in 



e 2 P - 1 

The numerical evaluation of the integral (we thank D. Volin for the help with that) is 
compatible with hi = h% = • • • = 0. Therefore 

Kp) = (B.6) 



which implies 

e 



:= / dxe^ px+lf{x) = -Z7T g2p _ i . (B.7) 



Applying the shift equation ( p.5|) we get (|B.4j) . 

Equation ( p.4|) is actually a particular case of the formula for the Fourier transforma- 



tion of the quantum dilogarithm presented in the Appendix of ref. |37j 
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